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ABSTRACT 



We compare the central mass distribution of galaxies simulated with three different models of the 
interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to I0 4 K, 
additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H2) with 
shielding of atomic and molecular hydrogen in addition to metal line cooling. The latter model 
includes a non-equilibrium calculation of the H2 abundance, self-shielding of H2, dust shielding of 
both HI and H 2 , and H 2 -based star formation efficiency. In order to analyze the effect of these 
models, we follow the evolution of four field galaxies with V pea k < 120 km/s to a redshift of zero using 
high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological ACDM context. 
CN| . Artificial observations were created for the galaxies and compared with real observations. The spiral 

galaxies produced in simulations with either primordial cooling or H 2 physics have bulge magnitudes 
and scale lengths very similar to observed galaxies and realistic, rising rotation curves. In contrast, 
the metal line cooling simulation produced galaxies with more massive and concentrated bulges and 
with the peaked rotation curves typical of most previous ACDM simulations of spiral galaxies. The 
Vh \ less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling 

or H 2 are linked to changes in the angular momentum distribution of the baryons. These galaxies 
had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar 
feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback- 
heated gas cools more slowly than when metal line cooling is included and so requires less energy to 
be expelled. When H 2 is included, the accompanying shielding produces large amounts of clumpy, 
cold gas where H 2 forms. Consequentially, star formation takes place in denser and clumpier gas. 
Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency 
of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material and 
formation of realistic bulges. The inclusion of the H 2 model and the resulting increase in feedback 
enabled the creation of simulations with both realistic metal line cooling of the halo gas and realistic 
structure within the central kiloparsec. 
Subject headings: 

1. INTRODUCTION ter (DM) halos was transformed into the stellar popu- 

The central question in A-Cold Dark Matter (ACDM) lations we observe "J galaxies today. _ To answer this 
galaxy formation is how the gas collected into dark mat- question, we must understand how the interplay between 
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the physical processes connected with star formation re- 
sulted in the observed stellar content of modern galaxies. 
These processes include the rate at which gas cools, how 
it is transformed into stars, and the regulation of star 
formation by such things as stellar feedback. Histori- 
cally, it has been very difficult for ACDM simulations 
to reproduced observed stellar bulges and the central re- 
gions of rotation curves. In this work, we concentrate on 
how different models of the interstellar media (ISM) and 
star formation impact the amount of mass in the central 
region of the disk galaxies, an area particularly sensi- 
tive to changes in the angular momentum distribution 
(AMD). We demonstrate how substantial redistribution 
of the angular momentum in simulated disk galaxies can 
be induced by changes to the ISM and, consequentially, 
to the star formation and supernova (SN) feedback. 

Simulations offer the ideal format for studying 
the evolution of the non-linear processes that domi- 
nate galaxy formation over cosmic time, such as the 
formation of galactic disks. Disks are created by pre- 
serving angular momentum from c osmic torques during 
the inf all of gas (IWhjtd fl98i [Barnes fc Efstathioul 
19871 iQuinn fc Binnevl Il992t Waller et al.l I2002t 
Vitvitska et al.l 120021) . If the gas infall is smooth, 
high- angular momentum gas will f orm a similarly 
high-angular moment um disk (e.g. iMo et all 119981 : 
iDalcanton et al.l I1997T) . In order to preserve the disk, 
large amounts of high-angular momentum material must 
be maintained throughout the galaxy's evolution. 

Artificial angular momentum loss has been shown 
through galaxy simulations to result from insuffi- 
cient feedback; gas cools and forms stars too quickly 
and angular momentum i s sub s equently lost dur- 
ing mergers tiWhite fc Reesl 119781; iDekel fc Silkl 119861: 
Navarro & Bcnz 1991; Efstathioul Il992t iKatz et al " 



1991 INavarro fc White! Il994 t iThoul fc Weinberg 



1997 



1996 : IQuinn et al l 119961 IN avarro fc Steinmeta 



ISommer-Larsen et al.l 119991: Gnedi 



Mailer fc Dekel 



ISommer-Larsen et aT 
2003- iRobertson et al 



van den Bosch et al.l 

Okamoto et al.1 120051: ID'Onghia et al.1 120061: iHoeft et al 




2006 1 : IGovernato et al l " 1200- 
20081: iScannapieco et aTT |20pj| 



Okamoto et al. 
Zavala et al.l 1200? 



Ceverino fc Klypinl 20091: Dutton fc van den Bosc 



20091: IKeres et al.l I2009t ISales et al.l 12010ft ." 

work has shown that angular momentum 



Further 
loss is 

caused by low resolution through such means as two- 
body he ating and artificial vis cosity dThacker et al 



2001 



20001: ISommer-Larsen et all 119991: IMaver et al 



IGovernato et all 12004 iKaufmann et al.l 12007 



Governato et al.l 120071: iNaab et all 120071: IMaver et al I 
20081: IPiontek fc Steinmetzl 1201 Other simulations 
have suggested that disk survival could be aided by 
reducing star form a tion efficiency at high redshift 
(|Gnedin fc Kravts"ovj 12010 V Recently, simulations 
have progressed in resolution and the modeling of 
feedback to the point that angular momentum is con- 
served and disks of the appropriate size that lie along 
the Tully-Fisher relation are produced even with active 



2006: iScanna 


pieco et al.l 20081: [Hopkins et al.l 12009b; 


Agcrtz et al. 


2010t Brooks et al.l 120111: iBrook et al.l 


2011aHGuedes et al.ll2011l). 



momentum imparted by tidal torques can be main- 
tained, the resulting galaxies do not match the mor- 
phologies of observed galaxie s. Although observed 
disk sizes can b e repr o duced (|Fall fc Efstat hioul 119801: 
Dalcanton et al.l 119971: IMo et al.l 119981: IBrooks et al.l 
20111 ). the resulting bulges predicted in A CDM are much 
too l a rge compared to observed bulges (|Bullock et al.l 



20011: Ivan den Bosch et~aH 120011: Ivan den Bosch! 12001 



Binney et al.l 120011; IMaller fe DekeH l2002t IMaller et al l 
2002t Ivan den Bosch et all 120021: ID'Onghia et al.l 120061: 
D'Onghia fc Navarro! 120071 ) . This disagreement is true 
across all disk galaxy masses, not only at the bulgeless 
dwarf galaxy scale at which the problem is most severe 
(|Duttonl [2009) . The discrepancy in the amount of low- 
angular momentum material argues that some process 
must alter the AMD of disk galaxies by preferentially re- 
move the low-angular momentum gas that would natu- 
rally builds up the central regions and bulges of galaxies. 

The most sensitive area in the disk to the amount of 
low-angular momentum material is the central 1 kpc. In 
this work, therefore, we focus on how the problem of 
overly large bulges in simulations is affected by different 
gas cooling models. The study of the central regions of 
galaxies in relation to angular momentum has only re- 
cently become attainable because of computational de- 
mands. Cosmological simulations are necessary to fol- 
low the detailed evolution of the AMD. Such simulations, 
however, are very computationally expensive as they re- 
quire force resolutions small enough to resolve the inner 
^lkpc, i.e. < 200 pc so that the inner lkpc contains at 
least five resolution units. 

A few such simulations that resolve the internal struc- 
ture of disks have demonstrated recently that low- 
angular momentu m gas can be preferent i ally removed 
from t he galaxy (IGovernato et all 120101: IBrook et aT 
l20T7al l20ia IMaccio et all l2012jT ~ 
120101) and |J 



1 120121: IMaccio et all 120121) . In IGovernato et al 
(J and IBrook et all (|2011al) . these outflows led to 



Importantly, though, even if the initial angular 



bulgeless disk g alaxies with realisti c rotation curves 
(|Oh et al.ll20Tll ). iGuedes et all (|201lD outlined the cre- 
ation of a Milky Way-mass galaxy with a realistic bulge- 
to-total disk mass r atio, created by centralized gas loss. 
IBrook et al.l (|2011bD demonstrated that the outflows cre- 
ated in these high resolution simulations naturally drive 
low-angular momentum gas from the galaxy. 

The preferential removal of low-angular momentum 
material is key to correctly simulating the central mass 
distributions in galaxies. In order for this removal to 
take place, it is crucial to realize the 3-D structure of 
the ISM at the scale of star formation. When star for- 
mation in simulations is limited to regions with densities 
comparable to those of observed molecular clouds, gas 
is able to collect in clumps with densities >100 amu/cc. 
When the massive stars in these clumps produce SN, 
the feedback is highly concentrated, leading to expand- 
ing bubbles of hot gas that puncture the disk and natu- 
rally drive galactic outflows. Importantly, the ability to 
resolve these dense star forming clumps allows the gas 
to reach comparable densities to the DM in the central 
regions of galaxies. Due to these comparable densities, 
sudden gas loss caused by SN feedback can flatten the 
gravitational potential, leading to an expansion of the 
DM orbits and a flattening of the DM density profile (i.e., 
a shift from a more cuspy p rofile to a more cored profile 
iPontzen fc Governatol l20lll) . Hence, both the creation 
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of outflows that remove low- angular momentum gas and 
the creation of DM cores rely strongly on star formation 
happening in regions of high density gas. 

A natural and physically-motivated method for 
ensuring that star formation and feedback take 
place in high density environments in the simula- 
tions is to connect star formation to the radiation- 
shielded areas of the ISM where H2 is abundant 
(|Robertson fc Kravtsovl 120081 ). In observed galax- 
ies, star formation is strongly co r related with molec- 
ular gas (iRownd fc Yound 11999} iMurgia et all l2002t 
Hever et al.ll2004l;|Gao fc Solomodl2004tlKennicutt et al l 
20071 [Blanc et al.ll2009t: iWarren et alj|2qi0t iBigiel et al.l 
20ld ISchruba et all 1201 U iBolatto et alj 120111) . Basing 
star formation on the H2 abundance also introduces ad- 
ditional dependencies, such as on metallicity. Because of 
this metallicity dependency, connecting star formation 
to H2 has been shown to lower the star formation effi- 
ciency in simulations of the low-metallicity , high-redshift 
Unive rse dGnedin fc Kravtsovl [2010l I20TI iKuHen et all 
I2012t iKrumholz fc Dekelll2012D . Reducing high-redshift 
star formation may lower t he amount of angular momen- 
tum loss during mergers ([Hopkins et al.l l2009bl ). How- 
ever, reducing early star formation cannot be responsi- 
ble for reducing the amount of low- angular momentum 
material below the initial AMD of the gas. 

As of yet, most simulations that have incorporated 
H2 based star formation have only used inefficient ther- 
mal SN feedback, making them unable to drive outflows 
that redistribute t h e ang ular momentum of the disk gas. 
iChristensen et all (|2012l) . however, demonstrated that 
using both H2-based star formation and efficient SN feed- 
back to simulate a dwarf galaxy changed the distribution 
of star formation and the temperature and density dis- 
tribution of the gas. In this paper, the addition of shield- 
ing with the H2 physics was largely responsible for the 
increased amount of dense, cold gas. Furthermore, an 
H2-based star formation recipe limited star formation 
to such dense gas. This ability of H2 to influence the 
location and environment of star formation suggests it 
could also affect the redistribution of angular momentum 
through feedback, making it a promising tool to better 
match observational trends in galaxies. 

The possible effect of the inclusion of shielded H2 
gas on feedback efficiency reflects a more general de- 
pendency of star formation and the feedback efficiency 
on the ISM model. The rates of cooling and heating 
naturally affect the mass, distribution, and dumpiness 
of the gas disk in galaxies. The ability of gas to ac- 
crete onto the disk of the galaxy regulates the amount 
of gas in the disk and, conseq uentially, the amount of 
gas available for star formation flChoi fc Nagaminell2009l ; 
ISchave et aLll201CH: Ivan de Voort et al.H201lD . Additional 
low-temperature cooling in the disk leads to the forma- 
tion of more clumps, increasing the amount of outflowing 
gas and further altering the AMD. On the other hand, 
greater amounts of cooling can als o make it more diffi - 
cult to expel gas through fe edback. IGuedes et all (|2011D 
and lGovernato et all (|2010f ) have produced galaxies with 
realistic disk scale lengths. However, both these galax- 
ies were simulated without full metal line cooling, rais- 
ing the question of how they would fare with a more 
realistic cooling mod e l. | Scannapieco etall ((20081 ) and 
iPiontek fc Steinmetzl ([201 ll ) allowed for cooling below 



10 4 K from metal lines but produce d galaxies wit h un- 
physically large bul ges. iBrook et all ([2011aL I2012I ). and 
iMaccio et all ([2012D also included metal line cooling in 
their simulations but were able to produce realistic bulges 
only through the addition of an extra source of energy 
from young stars. 

Here, we analyze the effect of different ISM models 
on star formation and the AMD in simulations of galax- 
ies. To do this, we simulated a small set of dwarf and 
spiral galaxies of different masses to z = using three 
different models of the ISM: primordial H+He cooling 
down only to 10 4 K, the addition of metal line cooling, 
and both metal line cooling and cooling from shielded 
H2 gas. Covering multiple galaxy masses and different 
observable quantities is necessary to constrain the mul- 
tivariate problem of the interaction between the ISM, 
star formation and feedback. We expect that the differ- 
ent ISM models will have less of an effect on the dwarf 
galaxies because of their smaller metallicities and lower 
H2 content. The inclusion of the additional cooling from 
metals or shielding will, therefore, have less effect on the 
amount and distribution of cold gas. As the galaxy mass 
increases, though, they will become more metal rich with 
larger amounts of cold, shielded gas. The choice of ISM 
model will, consequentially, have a greater effect on the 
structure of the ISM and the resulting efficiency of SN 
feedback. In order to compare our simulations to obser- 
vations, we generated and analyzed both optical and IR 
images as well as HI velocity cubes. We then examined 
several different metrics, including the rotation curves, 
and the characteristics of the bulges, the structure of the 
ISM, and the efficiency of feedback to determine the ef- 
fect of the ISM models. 

The paper is organized as follow: §2] discusses the sim- 
ulations; §3.11 describes the global observational prop- 
erties of the sample; §3.21 shows the differences in the 
mass distribution resulting from the introduction of H2 
related star formation; H3 31 describes differences in the 
ISM caused by the introduction of H2; and §33] compares 
the feedback efficiency. The results are then discussed in 
§3 and §5] is the conclusion. 

2. METHODS: DESCRIPTION OF THE SIMULATION 

The simulations described here are part of a new set 
of high-resolution simulations aimed at studying the for- 
mation of field galaxies. Our sample of halos includes 
two dwarf galaxies (virial masses of a few 10 10 M Q ) and 
two spiral galaxies (virial masses of several times 10 11 
M Q ). The properties of the simulated galaxies are listed 
in Tabic 1. 

The s imulations were co mputed in a WMAP3 cos- 
mology (ISpergel et al.ll2007h : O =0.24, A=0.76, h=0.73, 
08=0.77. Our halo sample was selected from lower res- 
olution cosmological volumes and then re-simulated at 
much highe r resolution using the 'zoom-in' technique 
(|Katzlll992l ) . Adopting the zoom-in approach enabled us 
to have a significant number of high-resolution particles 
while still following the surrounding large-scale structure. 
There were more than one million DM particles within 
the virial radius of each galaxy at z=0. For the dwarf 
(spiral) galaxies the DM, gas, and star particle masses (at 
the time of formation) were respectively: 1.6 (13)xl0 4 , 
3.3 (27.0) xlO 3 , and 1.0 (8.O)xlO 3 M . The force spline 
softening was 87 pc (170 pc) for the dwarf (spiral) galax- 
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Fig. 1. — Gas surface density and simulated observations of a spiral galaxy (h986) simulated with metal line cooling and H2 physics. 
Each frame is 24 kpc across. Left: HI surface density (gray scale) overlaid with red contours indicating the H2 surface density. The red 
contours are shown for the following levels: log S h 2 = -2, 0, and 2 amu/cc. Middle and right: SDSS multicolor (r, g, and i) images of the 
same galaxy at z=0 generated with SUNRISE in which the galaxy is oriented face on and edge on, respectively. This galaxy has a reddened 
H-band bulge half light radius, r e , of 0.83 kpc and a photometric bulge-to-total ratio of B/T = 0.43. 



ies. The minimum smoothing length for gas particles was 
0.1 times the force softening. This resolution was suffi- 
cient to follow the formation of star forming regions as 
small as 10 4 M Q . We further discuss resolution in relation 
to star formation in £12.21 

The simulations were performed with t he iV-body SPH 
code GASOLINE (jWadslev et al.ll2004h with a force ac- 
curacy criterion of 9 = 0.725, a time step accuracy of 
77=0.195 and a Courant condition of ?7c=0.4. The simu- 
lations i ncluded a cosmic UV field modeled following an 
updated lHaardt fc Madaul (|1996l ). which partially sup- 
presses the collapse of baryons into the smallest halos 
ijHoeft et al.1l2006HGovernato et al.ll2007ft . Our adopted 
star formati on and SN s c hemes have been described 
in det ail in iStinson et al.l (|2006l) and iGovernato et al.l 
(|2007f ). Briefly, the "blastwave" feedback scheme is im- 
plemented by releasing thermal energy from SN into gas 
surrounding young star particles. The heated gas parti- 
cles have their cooling shut off for a time equal to the 
momentum-conserving phase of the SN blastwave. This 
time is typically a few million years and is a function 
of the local density and temperature of the gas and the 
amount of energy injected and the cooling . 

As described in the following subsection, star forma- 
tion was limited to cold gas in dense regions, with cri- 
teria depending on the ISM model (see §2.2j) . Only two 
main free parameters in the SN feedback scheme (the 
star formation efficiency and the fraction of SN energy 
coupled to the ISM) were fixed to reproduce the prop- 
erties of present day g alaxies over a range of masses 
(jGovernato et al.l 12007). Without further adjustments, 
this scheme has been shown to rep roduce the relation 
between metallicity and stellar mass (jBrooks et al.ll2007T; 
iMaiolino et al.ll20"08l ) and the a bundance of Damped Ly- 
man a (DLA) systems at z=3 (jPontzen et al.ll2008D . 

2.1. Models of the ISM 

We computed the simulations using three different ISM 
models of increasing complexity and their correspond- 
ing star formation recipes. In these models, the differ- 
ent heating and cooling rates and the presence or lack 



of shielding resulted in varying temperature and density 
distributions. Because of the different structures of the 
ISM, interaction between the ISM, star formation, and 
feedback was also altered. Here we outline each of the 
models. 

All the considered models of the ISM include a cal- 
culation of the non-equilibrium ion abundances of H 
and He. The redshift-dependent UV background ra- 
diation was responsible for both gas heating and pho- 
toionization. Gas cooling was calcu lated based on col- 
lisional ionization (lAbel et al.lll997f). radiati ve recombi- 
nation (|Blackj 119811: iVerner fc ; Ferlandl Il996f ). photoion- 
ization, bremsstrahlung, and H and He line cooling (jCenl 
[1992] ). 

In our simplest ISM model, the heating and cooling of 
the gas is due only to the above reactions in the neutral 
and ionized H and He. For the purposes of this paper, we 
label this the "primordial" ISM model. Because of the 
lack of low-temperature coolants in this model, cooling is 
inefficient below 10 4 K. Gas is, therefore, unable to reach 
lower temperatures and the formation of high density 
regions is suppressed by the extra pres sure support. 

Additional cooling from metal lines (|Shen et al.ll2010() 
is included in our second ISM model (labeled "metals" 
throughout the paper). With the addition of metal line 
cooling, gas is able to cool more efficiently and to tem- 
peratures below 10 4 K. These metal line cooling rates 
were computed for different gas temperatures, densi- 
ties, metallicities, and amounts of UV background as- 
suming ionization equilibrium and optically thin g as us- 
ing CLOUDY (version 07.02: iFerland et al.l[T998f ). The 
metal line cooling rates were introduced to model the in- 
tergalactic medium. In this regime it is safe to assume 
that the gas is optically thin and that the cosmologi- 
cal UV background dominates the radiation field. In 
the ISM, however, this assumption can result in an un- 
derestimate of cooling rates. In the ISM, the lack of an 
interstellar radiation field beyo nd the UV background re - 
sults in substantially less CII ((Christens en et aT] [2012'). 
the most important coolant in spiral galaxies. Further- 
more, the lack of shielding from photoionizing radiation 
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prevents the additional formation of cold gas. 

Our most sophisticated ISM model (labeled as the 
"H2 " model) includes the no n-equilibrium formation of 
H2 (jChristensen et all I2012T ) in addition to metal line 
cooling. H 2 forms primarily on dust grains and is disso- 
ciated by Lyman- Werner radiation from young stars. It 
is, therefore, limited to dense regions of the ISM where 
the gas is shielded from Lyman- Werner radiation. Our 
model follows the hydrogen chemical network for each gas 
particle, assuming that the dust fraction is proportional 
to the metallicity and calculating the incident Lyman- 
Werner radiation based on the emission from nearby star 
particles. The model includes self-shielding and shield- 
ing by dust to protect H2 from photo-dissociation; dust 
shielding is also included when calculating the rates of 
HI photoionization and photoheating. The inclusion of 
shielded HI and H2 strongly decreases the amount of 
heating from photo-dissociation and ionization and re- 
sults in the formation of a cold (100 K) ISM. Some 
additional lo w-temperature cooling is provided through 
H 2 collision (|Gnedin fc Kravtsovl[20TTI) . While the more 
accurate calculation of CII and other low-temperature 
metal coolants would increase the cooling rates, the for- 
mation of cold gas is primarily dependent on the presence 
of shielding. This model for the ISM is, therefore, an im- 
portant step toward the accurate modeling of the cold, 
shielded phase of the ISM where star formation takes 
place. 

2.2. Star Formation 

As described in IStinson et al.l (|2006l ) , star formation 
in these simulations occurs stochastically in gas that is 
both sufficiently dense and cool. For a gas particle to be 
eligible to spawn a star particle in this model, it must 
have a temperature less than T ma x and a density greater 
than Pmin- For eligible gas particles, the probability, p, 
of a star particle forming is a function of the dynamical 

time, t/orm; 

p = ^i(l_ e -*At/<,o™) (1) 

where Ttig as is the mass of the gas particle, Tn S f ar is the 
initial mass of the potential star particle, and c* is a star 
forming efficiency factor. 

We used different star formation parameters (c*, T max , 
Pmin) based on the ISM model. For both the primordial 
and metal ISM models, c* = 0.1. We set T max = 10 4 K 
and pmin =10 (100) amu cm -3 for the spiral (dwarf) 
galaxies in order to select for the gas corresponding to the 
cold ISM in these simulations. This high density thresh- 
old for star formation was intended to limit star forma- 
tion to gas with densities similar to molecular clouds. 
The high value for T max reflects the inability of gas in 
the primordial ISM model or zero metallicity gas in the 
metal line cooling ISM model to cool below 10 4 K. Star 
formation is, therefore, allowed to occur in gas substan- 
tially hotter than the ~100K gas of molecular clouds. 

When using the H2 model, we incorporated the abun- 
dance of H2 into the star formation calculation as, 
star formation is observationally linked to H2 (e.g. 
iBigiel et al.l I2010D . We modeled this connection be- 
tween star formation and H2 by making the star for- 
mation efficiency a function of the local H2 abundance: 



c* = 0.1A"h 2 , where Ih 2 is the fraction of hydrogen 
in the form of H2. By making the star formation effi- 
ciency a function of the H2 abundance, we eliminated 
the need for a high density threshold to restrict star for- 
mation to molecular clouds. We, therefore, set the den- 
sity threshold to the arbitrarily low value of p m in = 0.1 
amu. We further limited star formation in this recipe 
to gas that has undergone molecular cooling by setting 
Tmax = 1000 K. We tested the sensitivity of star forma- 
tion to both p m i n and T max when using H2 and found 
that changes of a factor of five to either threshold did 
not significantly alter the star formation. With this star 
formation recipe, we naturally limit star formation to re- 
gions that correspond best to those in observed galaxies. 
Both the inclusion of H2 and the H2-based star forma- 
tion recipe were te s ted on a cosmological dwarf galaxy in 
IChristensen et akl (|2012t ) and found to produce a bluer 
dwarf galaxy at z = with more extended star formation 
and clumpier gas. 

The high density thresholds for star formation in com- 
bination with the low temperature cooling allow gas par- 
ticles to become Jeans unstable. While this can occur 
with metal line cooling, the formation of cold gas when 
shielding is included mean that it is f a r mor e prevalent in 
the H 2 ISM model. iBate fc Burkertl (|1997l ) showed that 
when the Jeans length or mass is unresolved artificial 
fragmentation can occur. At this point star formation 
becomes dependent on the specifics of the star forma- 
tion implementation. Many simulators have addressed 
this by requiring that Jeans instability never occur, ei- 
ther by not including low temperature cooling or adding 
a pressure floor. Still others made Jeans instability itself 
the criteria for star formation. We allow Jeans instabil- 
ity to occur with the understanding that gas particles 
that are Jeans unstable will soon form stars and that 
the highly efficient SN feedback will heat the surround- 
ing gas, preventing further collapse. While this method 
does make the star formation dependent on the imple- 
mentation, Jeans instability is a natural occurrence and 
it would be equally artificial to prevent it. We see no 
evidence of artificial collapse in our simulations. We cau- 
tion, however, that these simulations should not be used 
to study phenomena occurring below our resolution limit, 
such as the mass-spectrum of molecular clouds. 

2.3. Mock Observations 0/ the Simulated Galaxies 

In order to properly compare the outputs from our 
simulations to real galaxies, we calculated the observ- 
able properties of the simulated galaxies. These proper- 
ties include the photometric images, photometric surface 
brightness profiles, star formation rate (SFR) indicators, 
and gas surface densities of the galaxies. The observed 
properties of the four main halos are listed in Table 1 and 
the methods for calculating them are described below. 

The central gala xies in each simulation were first iden- 
tified using AHF f|Gill et al.ll200l iKnollmann fc Knebd 
[200l AMIGA Halo Finder,). AH F adop t s ove rdensi- 
ties as a function of redshift from iGrossI ()1997t ). with 
an overdensity of ^100 with respect to the critical den- 
sity, Pcrit, at z —0. The environment density at z=0 
is typical of field halos with moderate over-underdensity 
Sp/ p = -0.3-0.3 measured over a sphere of radius 3 h _1 
Mpc. The virial mass, gas mass, cold gas mass (defined 
as the atomic and molecular hydrogen mass multiplied 
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by a factor of 1.4 to account for helium), and stellar 
mass of the central galaxies (defined as the stellar mass 
within the virial radius that is not within satellites) are 
listed in columns 1-4 of Table 1. The galaxies simu- 
lat ed with the FT ISM model were previously analyzed 
in (jMunshi et alJl2012l ). in which they were shown to lie 
along the observed stellar mass-halo mass relation. The 
asymptotic circular velocities of the galaxies are listed 
column 5 of Table 1 and the method used to calculate 
them is described in §3.11 

A direct comparison with observations is vital to judge 
the realism of the simulated galaxies. We used the 
Monte Carlo radiation transfer code, SUNRISE (v 3.6 
lJonssonll2006D to generate artificial photometric images 
(see Figure Q} and spectral energy distributions (SEDs) 
of the simulated galaxies. The magnitude and photomet- 
ric im ages of our galaxies in Sloan D igital Sk y Survey 
(SPSS lAdeTman-McCarthv et al.l[2006h bands, QolmiQi] 
(11961 bands, and Two Micron All Sky Survey (2MASS 
iCohen et a l. 2003) bands were calculated by using SUN- 
RISE to convolve the generated SED with the filter trans- 
mission curves. The SDSS i-band magnitude of each 
galaxy is listed in column 6 of Table 1 and the SDSS 
g - r color is listed in column 7 of Table 1 . 

We used the SUNRISE-generated photometric images 
to calculated the surface brightness profiles of the spi- 
ral galaxies and decompose them into bulge and disk 
components. In order to measure this bulge-disk decom- 
positions of the simulated spiral galaxi es, we executed 
the sa me analysis discussed in detail in iFisher fc Droryl 
(2008). We opted to use the face-on 2MASS H-band 
images, as infrared data are less sensitive to the ob- 
scuring effects of dust and the face-on orientation pro- 
vides the most accurate fit. For each simulated galaxy, 
we determined the surface brightness profile through el- 
lipse fi tting. Ellipse fitting was d one using the rou- 
tine of iBender fc Moellenhoffl (|1987t ) using isophotal el- 
lipses of 256 equally spaced points increasing in radius. 
The software returns surface brightness, isophotal center 
(in X & Y coordinates), major and minor axis length, 
and position angle of the ellipse. The radial size of el- 
lipses were optimized to maintain a roughly constant 
signal-to-noise across the profile. For every galaxy, we 
also fit isopho t al ell ipses to the image with the soft- 
ware of iLaueil (Il98l . which is a Fourier based ellipse 
fitting method. This fit was used as an independent 
che ck on the quality of th e isoph otal ellipse fitting from 
the IBender fc Moellenhofl (Il987fl rout ine. The iso photes 
from both the IBender fc Moellenhoffl (|1987f ) and lLauerl 
(1985) routines were averaged together and the result 
used for the bulge-disk decomposition. 

The bulge-disk decompositions were determined by fit- 
ting two components, a Sersic bulge plus outer exponen- 
tial disk, to the major axis surface brightness profile. The 
inner radius cut for the fit was chosen to be 200 pc (« 
the softening length for the spiral galaxies) , and the outer 
radius cut was set to the location where the profile broke 
from a smoothly varying expone ntial. Bars were dealt 
with using the same method as in IFisher fc DrOTvl (2008, 
2010;. in which isophotes that are affected by bars were 
isolated and masked from the fitting. The bulge half- 
light radius, bulge Sersic index, and the bulge-to-total 
light ratio are listed in columns 8-10 of Table 1. 

SUNRISE also enabled us to measure spatially- 



resolved SFRs of the galaxies from the 24/xm and FUV 
emission observational indicators. We conver ted these in- 
dicato rs into SFRs using the formula listed in lLerov et al.l 
(2008) in order to make an accurate comparison to the 
resolved Kennicutt-Schmidt rel ationship generated for a 
large set of nearby galaxies in iBigiel et aD (|2008t ) (see 
fc|3.3[) . In order enable direct comparisons with a larger 
sample of observed galaxies, we used Ha emission to de- 
termine the SFR over the entire disk of the galaxy (col- 
umn 11 of Table 1). 

Finally, we measured the gas surface densities of the 
main halos in a way co mparable with THIN GS (The HI 
Nearby Galaxies Survey. IWalter etaLll2008D observations 
through the use of simulated HI velocity cubes. We cre- 
ated a velocity cube for each galaxy by rotating it to 45° 
angle and convolving each gas particle with the smooth- 
ing kernel and HI fraction calculated in GASOLINE. In 
order to mimic typical THINGS resolution and sensitiv- 
ity, we proceeded as if the galaxies were at a distance of 
5 Mpc. We used 1.5" pixels and applied beam smearing 
with a full width at half max of 10" to the data. We 
divided the data into 128 velocity channels, which had 
widths of 2.6 km/s for the two spiral galaxies and widths 
of 1.3 km/s for the dwarf galaxies. Finally, we made a 
sensitivity cut and ignored all emission from cells below 
the 2cr noise limit for the THINGS observations, in which 
a = 0.65 milli-Jansky/beam. This velocity cube was 
used to create images of the HI surface density (Figure Q] 
and Figure [9]) and to compare to the resolved Kennicutt- 
Schmidt relation (£ !3.3j) . 

3. RESULTS 

3.1. Global Properties of the Galaxy Sample at z=0 

Here we examine how the different models of the ISM 
and their interaction with the SN feedback can change 
both the total mass and mass distribution in galaxies. 

The interaction between the cold ISM, star formation, 
and SN feedback can affect the final baryonic content of 
the galaxies at z=0 by changing the amount of gas ex- 
pelled by feedback or stripped from the galaxy. We ana- 
lyzed this effect by determining the location of the sim- 
ulated galaxies on the baryonic Tully-Fisher (BTF) rela- 
tion. This relation compares the total amount of baryons 
in the form of either stars or cold gas within a halo to the 
rotational velocity and is known to be observationally a 
very tight relationship. We comp ared the main ha los of 
our simulations to the BTF from iMcGaughl (|2005| ). fol- 
lowing their observational analysis as closely as possible 
(Figured]). We defined Vf to be the asymptotic circular 
velocity of the gas. In order to calculate Vf, we fit either 
an increasing or decreasing arctangent function to the ro- 
tation curve (not including radii smaller than twice the 
softening length). The circular velocities, in turn, were 
calculated from the mass enclosed at each radii. In all 
cases except one primordial ISM dwarf galaxy, the sim- 
ulated galaxies lie exactly on the observed values. How- 
ever, as we will demonstrate below, the metal ISM model 
leads to a larger central mass distribution and higher ro- 
tational velocity in spiral galaxies. As the total baryonic 
mass increased, Vf increased in such a way that the sim- 
ulations moved along the BTF relation. This constancy 
across parameters suggests that the BTF is not suffi- 
cient as a discriminator between the ISM models, and 
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Fig. 2. — Baryonic Tully-Fisher relationship for the main ha- 
los overplotted on observational data. Filled red circles represent 
simulations with H2 and metal lines cooling ("H2"); open black 
squares represent simulations with only metal line cooling ("met- 
als"); open grey diamonds represent simulations with primordial 
cooling. The blue crosses are observational data from IMcGaugbl 
(2005). Wc calculated the baryonic (stellar and cold gas) mass 
fo r our simu l ations in the following way, using the method used 
in IMcG augh (2005) for observed galaxies. The simulation stellar 
masses were calculated from Ms , using a mass to light ratio deter- 
mined by the B - V color. The cold gas mass was the total HI and 
H2 (when applicable) mass, scaled by a factor of 1.4 to account for 
helium. As in the observational analysis, we used the asymptotic 
circular velocity from the rotation curve for Vt. Our simulations, 
with the exception of one dwarf with the primordial ISM model, lie 
on the observed data. The tight fit indicates that all ISM models 
produce galaxies that have the same relationship between central 
mass and Vf as observed galaxies. 

additional diagnostics are necessary. 

In addition to comparing the total baryonic masses, we 
also compared the distributions of matter produced by 
the different ISM models. Figure [3] shows the circular 
velocity curves for two of the four main halos (a repre- 
sentative dwarf and spiral galaxy) simulated using the 
different ISM models. The other two galaxies follow the 
same trends in mass distribution. These curves show the 
circular velocity at each radii calculated from the total 
mass contained within that radii. 

We found that all of the ISM models produced contin- 
ually rising circular velocity curves for the dwarf galax- 
ies. The similarity of the rotation curves for the dwarf 
galaxies is expected given; the low metallicity and H2 
content of these galaxies would result in relatively small 
differences between the ISMs produced by the three mod- 
els. In contrast, the circular velocity curves of the spiral 
galaxies are very different. Only the primordial and H2 
models produced continually rising rotation curves. The 
two spiral galaxies simulated with the metal ISM model 
have a central peak. A comparison s ample of observed 
galaxies is available in Ide Blok et al.l (|2008l) , which pro- 
vides a set of Hi-determined rotation curves for galaxies 
from THINGS. The gentle rise of the circular velocity 



curves for the galaxies simulated with either the primor- 
dial or H2 model is comparable to rotation curves of 
the observed galaxies. The central peak of the galaxies 
with the metal model, however, indicates a much higher 
central concentration and a build-up of low-angular mo- 
mentum baryons in conflict with observed constraints. 

In addition to the rotation curves, the changes to the 
mass distribution are also apparent in the DM profiles 
(Figure [4]). The DM profiles generated with different 
ISM models are shown for a representative dwarf and 
spiral galaxy. The dwarf galaxies experienced only slight 
changes to their DM profiles and, in contrast to the spi- 
ral galaxies, the metal ISM model produced the small- 
est central concentrations in both h799 and h516. As 
with the rotation curves, the similarities between the 
DM profiles can be explained by the low metallicities 
and molecular content of these galaxies, which reduces 
the differences between the ISM models. Similarly to 
the rotation curves, the most dramatic changes to the 
DM profiles occurred in the spiral galaxies: the DM is 
most concentrated for the metal simulations and least in 
the H2 simulations with the primordial simulations ly- 
ing in-between. The relatively low concentrations of the 
H2 and primordial simulations compared to the metal 
simulation is consistent with the lack of peaked rotation 
curves. The changes to the DM profile indicate that the 
reduced circular velocities seen in Figure [3] are not only 
the result of reduced central baryonic mass but also to 
an alteration in the galaxies's DM distributions. 

The changes to the rotation curves and DM profiles 
denote changes to the AMD. Figure [5] shows the AMD 
for both the DM and baryonic mass in a representative 
dwarf and spiral galaxy. The baryons in the plot include 
all observable baryons, i.e., the stars and the HI and H2 
gas (if applicable) multiplied by a factor of 1.4 to account 
for the cold helium (helium species are tracked separately 
from hydrogen in the code) . The DM AMD was normal- 
ized by dividing by the maximum of the distribution, 
while the baryon AMD was divided by the maximum 
and then further scaled by the baryon fraction within 
the halo (so that a maximum value of 1.0 indicates that 
the halo has the cosmic baryon fraction). The use of 
the metal ISM model resulted in a strong peak at zero 
angular momentum in the spiral galaxy. This peak was 
caused by a relatively large amount of very low-angular 
momentum material. In addition to lowering the central 
peak, both the primordial and H 2 ISM models resulted 
in less baryonic material in the disks of the spiral galax- 
ies, indicating greater amount of baryonic mass lost from 
the simulations. As was the case for the rotation curves 
and DM profiles, the AMDs of the dwarf galaxies show 
much less variation between models. 

3.2. Photometric Decomposition of Galaxies at z=0 

The reduction of the total stellar mass and the central 
concentration of baryons in simulations with either the 
primordial or H2 ISM raises the question: how do the 
different models of the ISM affect the distribution of stars 
in the galaxies? To answer this question, we compared 
the photometric profiles of the spiral galaxies. 

Figure |5] shows the 2MASS H-band radially averaged 
photometric profiles generated by SUNRISE for a repre- 
sentative spiral galaxy simulated with the different ISM 
models. From this figure, it is apparent that the pri- 
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Fig. 3. — The circular velocity profiles of the galaxies. The left panel shows a representative dwarf galaxy (h799) and the right panel 
shows a representative spiral galaxy (h986). Red lines represent simulations with the H2 ISM model; black lines represent simulations 
with the metal ISM model; grey dashed lines represent simulations with the primordial ISM model. The vertical dotted line indicates the 
radius at twice the softening length. The circular velocity curves of the dwarf galaxies vary only slightly between the different ISM models. 
Simulations of the spiral galaxies with the H2 or primordial ISM models have rising circular velocity curves, in contrast to the steeply 
peaked circular velocity curves of the simulations with the metal ISM model. 




Fig. 4. — Dark matter density profiles of the galaxies. The left panel shows a representative dwarf galaxy (h799) and the right panel 
shows a representative spiral galaxy (h986). Red lines represent simulations with the H2 ISM model; black lines represent simulations with 
the metal ISM model; grey dashed lines represent simulations with the primordial ISM model. The vertical dotted line indicates the radius 
at twice the softening length. The central densities of the dwarf galaxies are similar across all ISM models. The central concentration of 
the spiral galaxies is lowest for the H2 ISM model, followed by the primordial ISM model, and highest for the metal ISM model. This 
pattern is similar to the changes to the circular velocities with ISM model. 



mordial and H2 models produce galaxies with smaller 
central surface brightnesses than the metal ISM model. 
The lower central magnitudes in each denote a less mas- 
sive bulge and a secondary effect is the expansion of r e 
larger. 

In order to quantify both the decreased magnitudes of 
the bulges and larger values of r e we decomposed the 
H-band profiles of both the spiral galaxies into an ex- 
ponential disk and Sersic bulge using the method out- 
line in ^2.31 The values of r e , the Sersic index, and the 
B/T ratio for the two spiral galaxies are listed in Ta- 
ble 1. We compare the bulge H-band magnitudes and 
r e for the spiral galaxies with the different ISM models 
to observed galaxies in Figure [7] The observational data 
shown is from Fisher et al. (2012 in prep). This observed 
galaxy sample includes spiral galaxies from Hubble types 



of Sa to Sd and, therefore, will cover the range of bulge 
properties. It covers a sim ilar range in B-V as la r ge sur - 
veys of galaxies, such as Ide Vaucouleurs et al.l (|1991|) , 
and preferentially samples bright and large disk galax- 
ies, in which bulges are more commonly found. Galax- 
ies simulated with the primordial and H2 ISM models 
have both dimmer and more extended bulges, as seen 
in the photometric profiles. The reduced magnitude 
of the bulges occurs in connection to a smaller mass 
of low- angular momentum baryons in these two simu- 
lations. This result is consis t ent w ith the observational 
work of IRomanowskv fc Falll ()2012[ ). which shows differ- 
ent amounts of angular momentum for galaxies of differ- 
ent Hubble types. 

The bulges produced with the primordial and H2 ISM 
models are more similar to the bulges of observed galax- 
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Fig. 5. — The AMD of the DM (dashed lines) and disk baryons (solid lines). Results for a representative dwarf galaxy (h799) appear 
in the left panel and results for a representative spiral galaxy (h986) appear in the right panel. Red lines represent simulations with the 
H2 ISM model; black lines represent simulations with the metal ISM model; grey lines represent simulations with the primordial ISM 
model. For both the dark matter and baryons, the distribution, P(s), is shown for the angular momentum, j, scaled by the mean angular 
momentum, (j). These distributions of scaled angular momentum are then divided by the maximum of the distribution. For the baryons, 
the P(s) is again divided by the cosmic baryon fraction to produce the final plot. The baryons used for calculating the angular momentum 
lie within the optical radii and consist of the stars and cold gas (HI and H2 gas scaled by a factor of 1.4 for helium). The changes to 
the AMD are consistent with the changes in the rotation curves. In the spiral galaxies, the use of either the primordial or H2 ISM model 
results in decreased baryonic mass in the disk and a smaller mass of low-angular momentum baryons compared to the galaxies simulated 
with the metal ISM model. In the galaxies with the metal ISM model, the AMD is slightly centrally peaked. The angular momentum 
profiles of the dwarf galaxies are much more similar, although the galaxy simulated with the primordial ISM models has slightly more disk 
baryons and the galaxy simulated with the H2 model has a slightly lower central concentration. 
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Fig. 6. — Radially binned, face-on photometric profiles in the 
2MASS H-band for a representative spiral galaxy (h986). The 
red line represents the simulation with the H2 ISM model; the 
black line represents the simulation with the metal ISM model; the 
grey dashed line represents the simulation with the primordial ISM 
model. Simulations with the primordial and H2 ISM models have 
reduced central magnitudes. 

ies than those produced by the metal ISM model. As 
can be seen from the observational data, the bulges with 
the largest r e are also more luminous (and hence more 
massive). Rather than move along this relation, varying 
the cooling models in the simulations moves the result- 
ing bulge properties perpendicular to this trend. Recall 
from iJ3.ll that models moved along the BTF relation, 
negating the power of the BTF relationship to discrim- 
inate between the different models. The perpendicular 
movement by these models relative to the observed bulge 
size-luminosity relationship give this relationship greater 



discriminating power. 

Bulges have older average stellar ages than disks, im- 
plying that a large fraction of their stella r mass was 
formed at comparatively e arly redshifts (|Allen et al.l 
l2006tlMacArthur et al.ll2009[) . Furthermore, delayed star 
formation results in mergers happening between more gas 
rich galaxies, which leads to the reduced dissipation of 
angula r momentum and smaller bulges (Hopkins et al. 
I2009a|) . Given these two effects, we might expect that 
simulates with smaller bulges would have had less high- 
redshift star formation. We, therefore, compare the star 
formation histories for the galaxies with different ISM 
models to see if they can account for the differences be- 
tween the bulges. 

Figure |S] shows the star formation history for one of 
the spiral galaxies. We find that both the H2 and pri- 
mordial ISM models reduce the total stellar mass of the 
spiral galaxies, as can also be seen by comparing the 
data in Table 1. However, the reduction of star forma- 
tion does not necessarily take place at earlier times. In 
both galaxies, the star formation in these two models is 
especially reduced compared to the metal ISM model be- 
tween 6 and 8 Gyrs. This decrease is likely the result of 
greater gas-loss from feedback in the primordial or H2 
ISM models. The reduced gas mass in the galaxy limits 
the amount of star formation in both of these models. 
The amount of gas loss is further discussed in i j3.4l 

The use of a ^-dependent star formation law 
has been put forth as a way to reduce early star 
formation by raising the effective density thresh- 
old for star formation in low-metallicity gas (e.g . 
Robertson fc Kravtsovl I2008t IGnedin fc Kravtsovl [2010: 
Krumholz fc D ckcl 2012). However, we do not see a sig- 
nificant reduction of the stellar mass prior to 3 Gyrs in 
the simulations with the H2 ISM model. After this point 
in time, the metallicities of the galaxies have increased 
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Fig. 7. — Half light radius (r e ) versus magnitude (Mh) of the 
bulge for the spiral galaxies plotted against observed data. The 
red filled circles represent simulations with the H2 ISM model. 
The black open squares represent simulations with the metal ISM 
model. The grey open diamonds represent simulations with the pri- 
mordial ISM model. Observational data is from Fisher et al, 2012 
(in prep) and marked by dark blue and light blue filled squares, 
representing observations of bulges and elliptical galaxies, respec- 
tively. The black error bars are for the observed data. Both the 
primordial and H2 models of the ISM produce dimmer and larger 
bulges, in closer alignment with the observed galaxies. The metal 
ISM recipe, in contrast, produces galaxies with an excessive amount 
of bulge stars. 
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Fig. 8. — Star formation histories of the representative spiral 
galaxy, h986, produced by each of the models. The red line repre- 
sents the simulation with the H2 ISM model; the black line repre- 
sents the simulation with the metal ISM model; the grey dashed 
line represents the simulation with the primordial ISM model. 
There is a strong decrease in the star formation between 6 and 
8 Gyrs for the primordial and H2 ISM models compared to the 
metal ISM model. The different models do not, however, have a 
strong impact on the amount of stars formed before 3 Gyrs. 



enough to r educe its effect of t h e star formation law. As 
discussed in iChristensen et al.l ((2012), star formation in 
our simulations is strongly regulated by feedback so the 
reduced star formation efficiency in low metallicity gas 
has a relatively small effect. 

3.3. Properties of the ISM 

In this section, we compare the properties of the ISMs 
produced with each of the models to understand how 
changes between them are connected to changes to the 
mass distribution of the galaxies. The gas surface den- 
sities for a spiral galaxy simulated with each of the ISM 
models are shown in FigureJS] In this figure, the logarith- 
mic HI surface densities at the THINGS resolution and 
sensitivity are shown in grey-scale. For the H2 model 
ISM, the logarithmic H2 surface density is shown in the 
red contours. This figure shows the variation in the gas 
distribution produced by the different ISM models. With 
increased amounts of cooling, the gas distribution be- 
comes both clumpier and more extended. 

First, we analyzed the star formation environments in 
each of the simulations. To do this, we computed the 
resolved K-S relation for the four galaxies produced with 
the three different ISM model s, and compared th em to 
the THINGS data outlined in iBigiel et al.l poToh (Fig- 
ure Ho]) . In finding the resolved K-S relation for our 
sample of galaxies, we closely followed the observational 
methods. The SFRs were calculated based on the 24/im 
and FUV emission from SUNRISE simulated observa- 
tions. The HI surface densities were taken from the ze- 
roth moment of the HI velocity cube, which was com- 
puted at the same resolution and sensitivity of a stan- 
dard THINGS observation. The H2 surface densities 
were found directly from the simulation output (i.e. we 
did not calculate it from the CO emission as was done in 
observed sample). For each of the galaxies, we divided 
the SFRs and total hydrogen surface densities within the 
optical radius into 750 pc x 750 pc bins. The results are 
shown in Figure [TO] The SFRs and total hydrogen sur- 
face densities of the spiral and dwarf galaxies with each 
of the different ISM models are represented by the col- 
ored contours. The grey-scale filled contou r map rep- 
resents the observed data from IBigiel et~aLl (|2010f ) . In 
general, the simulated data follow the same trend as the 
observed data. Because of their low-mass, both the dwarf 
and low-mass spiral galaxies have surface densities lower 
than Ejj ~ 10 M Q pc -2 and lie on the Hi-dominated side 
of the K-S relation. The metal ISM model may have a 
slightly steeper slope than either the observed galaxies 
or the two other simulations. The difference, however, 
is slight and unlikely to have caused the differences in 
central stellar mass. We, therefore, conclude that dif- 
ferences in the K-S relation between the simulations are 
not responsible for changes to the bulge. The similar- 
ity between the results for the different ISM models also 
shows that the K-S relation is not a strong discriminator 
between our ISM models or star formation recipes in this 
regime. 

We next compared the densities and temperatures of 
the particles that form stars in the each of the models. 
Figure [11] shows the phase diagram for the gas parti- 
cles immediately prior to forming stars. The addition 
of shielding of both HI and H2 to the simulation with 
the H 2 ISM model resulted in a build-up of gas at w 100 
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Fig. 9. — Gas surface density of a spiral galaxy (h986) simulated with each of the ISM models: left (primordial), middle (metals), right 
(H2). The gray scale is the HI surface density at the THINGS resolution and sensitivity, ranging from log Sjji = -2 to 2 amu/cc. In the 
right panel, the red contours show the H2 surface density at log Eh 2 = -2, 0, and 2 amu/cc. The dumpiness, amount, and extent of disk 
gas increase from left to right. 




Fig. 10. — Resolved K-S relation for the three different models of the ISM: left (primordial), middle (metals), right (H2). These plots 
show the SFRs and total gas (b oth HI an d H2) surface densities averaged over 750 pc X 750 pc bins. The grey-scale filled contours are 
observational data from Bigiel ct al. (2010) for local dwarf and spiral galaxies. The colored contours represent the combined data for the 
simulated spiral and dwarf galaxies. We calculated the SFRs from the SUNRISE-generated FUV and 24/^m emission, the HI emission from 
the face-on velocity cube of average THINGS resolution and sensitivity, and the H2 directly from the simulation. All the ISM models 
produced a star formation law on kpc scales that is consistent with the observed data, although the metals ISM model may have a slightly 
less accurate fit. 



K. Furthermore, the ^-dependent star formation recipe 
caused star formation to take place in denser gas with 
higher H2-to-HI fractions. These densities are frequently 
higher than the current density where gas transitions 
from HI to H2 because of the amount of high-redshift star 
formation that took place in lower metallicity gas and 
because the delay in star formation after H2 has formed 
allows for further collapse. In the metal ISM model, the 
addition of metal line cooling resulted in star formation 
taking place at temperatures lower than 10 4 K. However, 
gas in this simulation generally formed stars soon af- 
ter crossing the 10 amu cm~ 3 density threshold without 
cooling significantly first. Without additional sources of 
cooling, gas in the primordial ISM model formed stars at 
10 4 K. This comparison illustrates how the star forming 



gas may be similar when averaged over ~ 1 kpc scales 
(as in the resolved K-S law shown in Figure [T0|) but very 
different on the smaller scales. These changes to the den- 
sity of the star forming gas have the potential to affect 
the efficiency of the blastwave feedback model, with more 
dense gas resulting in a more concentrated energy. 

Changes to the ISM model affect not only the den- 
sities of the star forming gas but also the densities of 
nearby gas particles. The effect of the different ISM 
models on the gas density and azimuthally-averaged sur- 
face density across the disk of the galaxy are shown in 
Figure [T^l In this plot, the three panels represent the 
different ISM models. Within the optical radii there is 
little change in the average surface density across sim- 
ulations. However, the densities of individual particles 




Fig. 11. — Characteristics of the star forming gas in a repre- 
sentative spiral galaxy (h986) simulated with each of the different 
ISM models and star formation recipes. The contours show the 
logarithmic distribution of gas particles that formed stars in den- 
sity and temperature; they represent the number of star particles 
formed and are spaced linearly. The the red contours are for the 
H2 ISM model, the black contours are for the metal ISM model, 
and grey contours are for the primordial ISM model. Compared to 
the other two models, star formation in the H2 ISM model takes 
place in colder gas because of the addition of dust shielding of HI 
and H2 and denser gas because of the requirement for gas to be 
molecular to form stars. 
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Fig. 12. — Density and azimuthally averaged surface density of 
gas as a function of radii for a representative spiral galaxy (h986) 
simulated with the different ISM models The top panel shows the 
primordial ISM model, the middle panel shows the metal ISM 
model, the bottom panel shows the H2 model. Grey points cor- 
respond to the left axis and represent the density of gas particles. 
The black curves correspond to the right axis and represent the 
radially-binned average surface densities of the gas. Vertical lines 
mark the optical radii. The increased cooling from metals and the 
shielding in the H2 model resulted in clumpier gas and greater 
amounts of high density gas although the average surface densities 
were similar for the different simulations. 



Fig. 13. — Power spectrum of the HI surface density in a repre- 
sentative spiral galaxy (h986). Filled red circles represent simula- 
tions with the H2 ISM model; black empty triangles represent simu- 
lations with the metal ISM model; grey empty diamonds represent 
simulations with the primordial ISM model. The 2-dimensional 
Fourier transform was calculated on the face-on HI surface density 
and a power-law was fit to the power spectrum (solid lines). The 
introduction of H2 produces a shallower slope than the other two 
ISM models, indicating that the gas has more power at small scales 
(i.e. clumpier). 

vary substantially between the models. The increased 
cooling in the metal line model and the shielding in the 
H 2 model resulted in the gas particles reaching progres- 
sively higher densities. High-density star forming regions 
appear in the plot as spikes in the grey points where gas 
particles at similar radii hav e especially high densities. 
IPontzen fc Governatol ([201 ID showed that dense clumps 
of gas, such as these, are an important factor in the ef- 
ficiency of SN feedback at transferring angular momen- 
tum. When SNe explode, the energy from the blastwave 
feedback model is concentrated in a smaller number par- 
ticles than if the star forming gas (and the surrounding 
gas) was more diffuse. This concentrated burst of energy 
results in the affected gas particle being heated to higher 
temperatures and for longer. 

The formation of these clumps also allows for star for- 
mation throughout a greater extent of the disk. The 
slight increase in the optical radii (marked by the verti- 
cal lines) in the ISM models with more cold gas indicates 
a greater spatial extent of star formation. As shown in 
([Christensen et aT]|2012|) . the additional cold gas allows 
star forming regions to form farther out in the disk of 
the galaxy in the H 2 simulation. 

We quantified the dumpiness of the gas for each of 
the ISM models using the power spectrum of the HI 
surface density distribution (Figure H2J) ■ We calculated 
the power spectrum for the each of the main halos with 
the different ISM models by computing a 2-dimensional 
Fourier transform of the face-on HI surface density using 
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a Fast Fourier Transform algorithm. For the HI surface 
density, we used the first moment of the HI velocity cube. 
However, we did not apply the THINGS beam smoothing 
or make a sensitivity cut for this analysis as most obser- 
vations of the power spectra of external galaxies are done 
with resolutions comparable to or higher than our soft- 
ening length. We truncated the power spectrum at the 
force softening length (170 pc in the spiral galaxies). In 
all simulations, the power spectra can be fit with power 
laws. The power indicates the relative amount of energy 
at different scales with shallower slopes indicating greater 
the energy on small scales, i.e. greater dumpiness. The 
exponents of the power law fit for h986 (h603) are —2.8 
(—3) in the primordial cooling model, —2.6 (—2.7) in 
the metal cooling model, and —1.9 (—2.1) in the H2 
model. The fits for observed galaxies on scales greater 
than 200 pc are consistent wi th these values and have 
range between —1-5 and —3 (|Stanimirovic et al.l 119991: 
Elmegreen et al.l 120011: iBegum et al.l 120061: iDutta et al.l 
2008L 120091) ." Therefore, we cannot use the power spec- 
trum to distinguish between the accuracy of these mod- 
els. However, it is clear that the shielding included in 
the H2 model produced a clumpier gas with more power 
at small-scales. We see no strong differences between the 
power spectra of galaxies with the pr imordial or metal 
ISM models. iPilkington et al.l (|201lD . in contrast, de- 
tected substantially more power at large scales in a dwarf 
galaxy run with metal line cooling compared to a dwarf 
galaxy without it. That galaxy, however, was caught in 
a state without any on-going star formation, which may 
account for the enhanced formation of larger-scale struc- 
ture, such as spiral arms. 

3.4. Feedback- Induced Gas Loss 

These changes to the ISM have the potential to affect 
the distribution of gas throughout the galaxy by changing 
the efficiency of gas loss from feedback. In order to study 
the efficiency of feedback, we tracked gas particles over 
the history of the simulation. We followed the history 
of the most massive progenitor galaxy and traced when 
gas particles were accreted to it, formed stars, and were 
ejected. In order to identify the gas particles ejected from 
the galaxy by SNe, we first identified the gas particles 
that had been part of the disk of the galaxy (p > 0.1 
amu cm~ 3 , T < 2 x 10 4 K and with a vertical distance 
to the central plane of the galaxy of less than 10 kpc) at 
any output. We then looked for subsequent outputs in 
which the gas particle 1) had been heated by a SN since 
the time it was first identified as being in the disk and 2) 
was no longer part of the disk. These two criteria defined 
the gas particles ejected from the disk. We also defined 
a subset of the ejected gas to be the gas expelled from 
the main halo by selecting particles that traveled beyond 
the virial radius after being heated by SNe. This gas will 
only ever again become a part of the galaxy is the virial 
mass is substantially increased following the expulsion. 

This method of identifying outflows is limited by the 
timcstep resolution. The outputs for our simulations are 
spaced at approximately 350 Myr intervals. Therefore, 
this method does not identify particles that were ejected 
from the disk by SNe and then re-accreted onto the disk 
between two timesteps. Since the simulations all share 
similar time resolution, however, the results may be used 
to compare them, even if the accuracy of the absolute 
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Fig. 14. — Cumulative history of gas loss over the history of the 
representative spiral galaxy (h986) divided by the amount of gas 
ever part of the galaxy. Red lines represent simulations with the H2 
ISM model; black lines represent simulations with the metal ISM 
model; grey lines represent simulations with the primordial ISM 
model. Solid lines show the fraction of gas ejected from the galaxy 
disk and dashed lines show the fraction of gas expelled beyond 
the virial radius. The amount of gas ejected is defined to be the 
amount of gas what was once in the disk of the galaxy but in some 
subsequent timestep had been heated by SNe and was no longer 
within the disk. The expelled gas is the subset of the eject gas that 
traveled beyond the virial radius of the halo. In order to account 
for differing amounts of gas accretion, we divided the amount of gas 
lost by the total mass of gas particles ever contained in the galaxy. 
The greatest amount of mass loss happens in the primordial ISM 
galaxy, followed by the H2 galaxy, and, finally, the galaxy with 
metal line cooling. 

amounts of ejected gas is limited. Within our definition 
we do not identify the multiple times a gas particle may 
be ejected from the disk as separate instances but only 
count the first time. 

The mass of gas lost is shown in Figure [14] In this 
figure, the cumulative histograms of gas ejected from 
the disk and expelled from the main halo over time are 
shown for a representative spiral galaxy simulated with 
the three different ISM models. In order to control for 
slightly different amounts of gas accretion, the distribu- 
tions were normalized by the total mass of gas particles 
ever contained within the virial radius of the main halo. 
The results for the other spiral galaxy are comparable 
with the plot shown. The total fraction of gas lost from 
each of the four galaxies is listed in Table 1. 

The effect of the ISM on the spiral galaxies follows 
the same pattern as the AMD suggests: the the primor- 
dial ISM galaxy experienced the most gas loss from feed- 
back, followed by the H2 galaxy, and finally the galaxy 
with the metal ISM model. The spiral galaxies simulated 
with either the primordial or H2 models also have sub- 
stantially reduced central concentrations and substan- 
tially greater mass loss than the metal ISM model. This 
pattern supports the theory that feedback is responsi- 
ble for changing the central mass distribution of galax- 
ies. The difference between the "ejected" and "expelled" 
gas represents the fraction of gas particles that either re- 
main in the halo indefinitely or are re- accreted onto the 
galaxy. In the spiral galaxy, h986 (h603), 53 (50)% of 
the ejected gas in the primordial ISM model is entirely 
expelled from the galaxy, as opposed to 42 (38)% in the 
metal ISM model, and 36 (31)% in the H 2 model. The 
difference in the fraction expelled is caused by the ad- 
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Fig. 15. — Mass loading factor for a representative spiral galaxy 
(h986) simulated with the three different ISM models. The red 
line represents simulations with the H2 ISM model; the black line 
represents simulations with the metal ISM model; the grey line 
represents simulations with the primordial ISM model. The mass 
loading factor is defined in this plot to be the mass of gas ejected 
from the disk of the galaxy in 0.5 Gyr bins divided by the stellar 
mass formed during that same time period. 

dition of metal line cooling. This addition increases the 
rates of cooling for ejected gas particles in the halo, mak- 
ing them less likely to be completely expelled. Previous 
research with simulations has als o shown galactic winds 
can enrich the circumg alactic gas (|Oppenheimer Sz Pavel 
120061 : iFord et all l2012t ) and that a substantial amount 
of gas in galactic winds is r e-accreted at a later tim e 
(|Qppenheimer fc DavdfeOOfl: lOppenheimer et al.l [2010h . 
Despite our different feedback models and methods for 
identifying outflows, our simulations appear to be gener- 
ally consistent with these previous results. 

The total amount of gas ejected from the galaxy disk is 
a function of the SFR. In order to compare the efficiencies 
of SNe at ejecting gas in these simulations, we divided 
the mass of gas ejected by the amount of stars formed 
to arrive at a mass loading factor. The mass loading 
factor in 500 Myr time bins as a function of time for the 
representative spiral galaxy with the three different ISM 
models is shown in Figure 1151 In this figure it is clear 
that both the H2 and primordial ISM models resulted in 
greater mass loading than the metal model. This trend 
is as expected from the comparison of the mass profiles 
and bulges. 

4. DISCUSSION 

The correspondence between gas loss and central mass 
concentration in galaxies illustrates how stellar feedback 
can reduce the amount of low-angular momentum mate- 
rial in galaxies. The differences between the ISM mod- 
els are likely the result of two different competing ef- 
fects. On the one hand, greater cooling results in more 
cold gas, which requires greater amounts of energy to 
be expelled from the galaxy. On the other hand, in- 
creased gas dumpi ness has been shown t o lead to more 
efficient feedback (IGovernato et all 120101 : IGuedes et al.l 
1201 It IGovernato et al.ll2012| ). When star formation and, 
therefore, SNe occur in high density peaks in the blast- 
wave FB model, the energy from SNe is concentrated in 
a smaller number of gas particles than if the gas and star 
formation were evenly distributed throughout the disk. 



As a result, these particles are ejected farther from the 
disk and have a larger initial density. Therefore, they 
have a greater effect on the gravitational potential of the 
galaxy. The accumulation of these sudden changes to 
the gravitational potential are what cause the reduction 
of the central DM concentration over time. 

The comparatively high efficiency of gas loss from 
galaxies with the primordial ISM model was the result 
of the lower cooling rates compared to the other two 
models. Gas heated by SNe in this model both starts 
at higher temperatures (> 10 4 K) and cools more slowly 
than when metal line cooling is included. The result of 
this decreased cooling in our simulations of spiral galax- 
ies was a greater amount of gas ejected farther from the 
disks. When metal line cooling is added, the tempera- 
ture of the star forming gas decreases and the cooling 
rates increases. As a result, greater amounts of energy 
are required to eject the gas and once ejected, the gas 
cools and re-accretes onto the disk more quickly. In 
our simulations, these effects resulted in spiral galaxies 
with more low-angular momentum material and higher 
central concentration. A smaller fraction of the gas re- 
mained outside of the disk long enough to be identified 
and being ejected and a smaller fraction was expelled 
from the galaxy, resulting in smaller mass loading frac- 
tions. Other simulations with metal line cooling have 
anticipated the effect of the increased cooling on feed- 
back and compensated for it b y increasing the amoun t 
of energy released per SN (e.g. Pilkington etall \20TS} . 
A full comparison such as ours, in which the energy re- 
leased per SN is held constant, is necessary, however, to 
find the extent to which the additional cooling affects the 
matter distribution in simulations. 

When H2 is added, the inclusion of the dust shield- 
ing for both HI and H2 and self-shielding of H2 results 
in larger amounts of cold dense gas. In our simulations, 
the formation of this cold, shielded gas greatly increased 
the dumpiness of the ISM, as was evident when com- 
paring the gas densities of the particles and the power 
spectrum of the HI surface density. Furthermore, by 
linking star formation to the H2 abundance, star for- 
mation was forced to occur in much denser gas than in 
either the primordial or metal ISM models. Together, 
the shielding of gas and the ^-dependent star forma- 
tion recipe forced star formation to occur in very dense 
clumps of gas. While the temperature of the star form- 
ing gas was decreased and while the metal line cooling 
models were still included, the increased dumpiness still 
produced more efficient mass loss from SN feedback in 
the spiral galaxies. The concentration of SN energy into 
the high density, star forming clumps was sufficient in 
these galaxies to increase the mass loading factor for the 
ejected gas up to that of the primordial cooling model. 
Despite the increased amount of ejected g smaller 
fraction of that ejected gas than in the primordial model 
was fully expelled from the galaxy because of the metal 
line cooling of halo gas. This increased feedback effi- 
ciency resulted in decreased central concentration and 
smaller bulges than with the metal ISM model. 

These results raise the question: what aspects of the 
H2 ISM model primarily caused the differences in the 
simulations? The first change was the increased amount 
of cold gas generated by the dust and self-shielding, 
which resulted in a clumpier ISM. The second change 
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was forcing of star formation to occur in shielded, dense 
gas, where the H2 abundances were high. To first order, 
both of these effects could be reproduced by implement- 
ing a dust shielding model for HI and a star formation 
recipe that was a function of the local shielding. Star for- 
mation laws that are a function of the shielded gas (as op- 
posed to the H 2 ) has been proposed (iMac Low fc Gloverl 
|20H IKrumhohd l2012t iGlover fc Clarkl 12012ft as a phvs- 
ical mechanism to explain the observed connection be- 
tween H2 and star formation. It is possibly that the time 
delay introduced by requiring shielded gas to form H2 
may result in subtle differences between H2-based and 
shielding-based star formation laws. In the future, de- 
tailed simulations may be able to tease the effects apart. 

In contrast to the spiral galaxies, the dwarf galaxies 
show much less variation in their profiles and gas loss. 
The smaller dependency of the dwarf galaxies on the 
ISM model is the result of their lower metallicity. Lower 
metallicities result both in less metal line cooling and less 
molecular hydrogen. Therefore, the resulting ISMs are 
more similar across the different models. The contrast 
between the dwarf galaxies and spiral galaxies illustrates 
the importance of testing galaxies of different masses. 

Our analysis also demonstrates how two of the more 
common metrics for analyzing galaxies, the Tully-Fisher 
relation and the K-S relation, are not sufficient for com- 
paring galaxy models. We found that our galaxies moved 
along the Tully-Fisher relation as both the mass of the 
disk and the circular velocity of the disk increased to- 
gether. The Tully-Fisher relation was, therefore, insensi- 
tive to significant changes to the total baryonic mass and 
central mass distribution. Both the rotation curves and 
the bulge properties were much more sensitive probes of 
changes to the AMD and provide stronger observational 
constraints. 

We also found that our very different ISM models and 
star formation criteria were all able to produce galaxies 
that all lay along the resolved K-S relation. Even when 
measured over 750 pc squares, small-scale variations in 
the ISM were washed out. Furthermore, in these simula- 
tions feedback is such a strong regulator of star formation 
that differences between the star formation thresholds 
and star formation efficiency had a relatively small ef- 
fect. Wc found the gas particle densities and HI surface 
density power spectrum to provide a more sensitive com- 
parison of the structure of the ISM produced with each 
of the models. 

5. CONCLUSIONS 

We analyzed a set of high-resolution dwarf and spi- 
ral galaxies computed with three different ISM models: 
1) primordial cooling only, 2) metal line cooling and 3) 
metal line cooling with gas shielding, non-equilibrium H2 
formation, and an H2-based star formation recipe. These 
galaxies were simulated in a cosmological context to a 
rcdshift of zero. In order to make our results most com- 
parable to observational data, we analyzed the stellar dis- 
tribution and mass and the density structure of the ISM 
at z = in an observationally motivated fashion through 
the creation of simulated optical images and HI velocity 
cubes. We determined the effect of different ISM mod- 
els on the efficiency of feedback in redistributing matter 
throughout the galaxy and we successfully created spiral 
galaxies with reduced amounts of low- angular momen- 



tum material. We found the following results for our 
suite of simulations: 

• Except for one dwarf galaxy simulated with the pri- 
mordial ISM model, all galaxies lay along the BTF 
relation. The changes to the ISM simply caused the 
galaxies to move along the BTF relation, making 
it poor discriminator between these models. 

• The different ISM models resulted in clear changes 
to the centers of the spiral galaxies. Spiral galaxies 
simulated with both the primordial and H2 ISM 
model had reduced masses within the central lkpc 
compared to the metal ISM model, which was ap- 
parent in their realistic non-centrally peaked ro- 
tation curves, less concentrated DM, and reduced 
mass of low-angular momentum material. 

• From fitting Sersic and exponential components to 
H-band images of our spiral galaxy simulations, we 
found that both the primordial and H2 models of 
the ISM produced realistic bulges, whereas the sim- 
ulations with only metal line cooling produced too 
massive of bulges. 

• The star formation histories of the spiral galaxies 
with the metal ISM model showed increased star 
formation around a redshift of 1 compared to the 
other two models. However, there is little difference 
between the star formation histories for the spiral 
galaxies prior to 3 Gyrs. This similarity at high 
redshift suggests that the smaller bulges produced 
with the H 2 ISM model were not the result of re- 
duced star formation efficiency in low metallicity 
gas. 

• All simulated galaxies fell along the observed K- 
S relation (resolved to bins of 750 2 pc 2 ). How- 
ever, the individual star forming particles had very 
different density and temperature distributions for 
the different models. In particular, in the H2 ISM 
model stars formed from colder and denser gas be- 
cause of the presence of shielding and the use of an 
H2-based star formation law. 

• Increasing amounts of cold gas resulted in clumpier 
ISMs. The effect was particularly strong for the 
spiral galaxies with H2. These galaxies had shal- 
lower power structures and increased amounts of 
high density gas. The ability of gas to collect in 
high-density clumps at larger radii resulted in more 
extended star formation. 

• The spiral galaxies simulated with the primordial 
and H 2 ISM models had greater amounts of gas 
ejected beyond the disk by SNe and larger mass 
loading factors. This increased efficiency of feed- 
back was likely responsible for the changes to the 
central kpc of the galaxies. The galaxies with the 
metal and H2 ISM models had smaller fractions 
of their ejected gas escape beyond the virial radius 
because of increased cooling in the halo from metal 
lines. 

• The dwarf galaxies simulated with different ISM 
models had relatively similar structures, includ- 
ing rotation curves and stellar profiles, because 
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of their smaller metallicities. In these galaxies, 
there was less potential for metal line cooling and 
shielding/ H 2 , which decreased the effect of the dif- 
ferent ISM models. 

The fact that both decreased cooling through lack of 
metal line cooling and increased cold gas from the shield- 
ing in our H2-model resulted in decreased bulges and 
increased feedback efficiency implies that multiple pro- 
cesses are at work in determining the mass of bulge stars 
formed. Compared to galaxies formed with only pri- 
mordial cooling, the additional cooling in the metal ISM 
model meant that more feedback energy was required to 
expel particles. The result was spiral galaxies with more- 
cuspy DM profiles, greater amounts of low-angular mo- 
mentum gas, centrally-peaked rotation curves, and big- 
ger and more concentrated bulges. With the addition 
of our H 2 model, the efficiency of feedback once again 
increased. In this model, star formed within extremely 
dense clumps, resulting in a greater concentration of SN 
energy and a greater likelihood of ejecting gas particles 
from the disk. The increased efficiency of feedback re- 
sulted in spiral galaxies with smaller central DM concen- 
trations and decreased amounts of low-angular momen- 
tum material, which led to more realistic rising rotation 
curves and smaller bulges. This final model is able to pro- 



duce realistic galaxies by including the metal line cooling 
necessary for modeling the CGM and the re-accretion of 
gas onto the disk and by including a model for shielded, 
molecular gas and the star formation that raised the feed- 
back efficiency enough to produce the observed centers of 
galaxies. 

Forming realistic centers of galaxies requires reducing 
the amount of low-angular momentum material in galax- 
ies. Together, these simulations reinforce the idea that 
stellar feedback is an effective way to remove low-angular 
momentum material from galaxies and to shape the cen- 
tral DM distribution. They also demonstrate that stellar 
feedback models are inexorably linked to ISM models in 
simulations. In order to create realistic galaxies, models 
of feedback, ISM and star formation must work together 
in a realistic manner. 
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Note. — Properties of the main halos at z —0. The virial mass, gas mass within the virial radius (R v i r ), cold gas mass, and stellar mass arc calculated directly from the simulations. 
The cold gas is defined as being 1.4(Mhi + A/h 2 )i where Mm is the total mass of atomic hydrogen, Mh 2 i s ^he total mass of molecular hydrogen, and the factor of 1.4 compensates for 
the mass of cold He (He species arc tracked separately from hydrogen in the code). In the primordial and metal simulations, the H2 mass is zero and the hydrogen gas that would be 
molecular remains as HI. Vf was calculated for each galaxy by fitting an arctangent function to the circular velocities. The SDSS i, 3, and r band magnitudes were calculated for a disk _ 
inclination of 45° from a simulated SUNRISE observation, which includes dust obscuration. SUNRISE was also used to determine the H-band profile. For the two spiral galaxies (h603 ^ 
and h986), a Sersic bulge and exponential disk were fit to the H-band profile. The bulge half light radius (r e ), the Sersic index of the bulge, and the B/T ratio for these fits are listed p. 
in columns 8 — 10. The two dwarf galaxies (h516 and h799), in contrast, were best fit with only a n exponential disk. The SFR at z — for each of the galaxies was calculated from the ^ 
SUNRISE-gcnerated H a emission, using the same conversion as lKennicutt Jr. Kennicuttj fll998T) . The "% Gas Lost" refers the the mass of supernova-heated gas expelled beyond the 2 
virial radius divided by the total mass of gas particles ever within the virial radius of the galaxy (see §3.4 for details). cc 
a Appears with primordial cooling and an additional low temperature cooling mode 1 in[Govci'iiato ct al. (2010). The galaxies arc labeled as DGl(h516) and DG2 (h799). P 
b Appears with the Primordial ISM mo del in [Brooks et all ^OlH ). In lBrooks et all i2011|) . both h603 and h986 arc at lower resolution. ^ 
c Appears with the Metal ISM model in I Oh et all J2011f) as DGl (DG2). p 
d Appears with the Metal ISM model in Pontzcn &c Govcrnato (2011). 
e Appears with the H 2 ISM model in Christcnscn ct al. (2012). 
^ Appears with the H 2 ISM model in[Govcrnato ct al. (2012). 
£ Appears with the H 2 ISM mode 1 in lMun shi ct al. [ 20121) . 



